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Wilson fermions quark bilinears to three loops 
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Quark currents renormalization constants can in principle be safely computed in lattice pertur- 
bation theory. In practice, traditional lattice perturbative computations are quite cumbersome, 
so that so far only the first loop results were available. By making use of Numerical Stochastic 
Perturbation Theory we reached three (and with less statistical precision even four) loops, both 
in quenched and in unquenched theory. Convergence properties of the series can be assessed and 
comparison with non perturbative results (where available) can be made: high loops computations 
of renormalization constants can be a valuable tool for lattice QCD. 
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1. Introduction 

Despite the fact that logarithmic divergent renormalization constants like those for quark bi- 
linears could be in principle safely computed in Lattice Perturbation Theory (LPT), one usually 
tries to compute them non-perturbatively, relying on some intermediate scheme. Popular choices 
are RI'-MOM [jl|] and SF Drawbacks of traditional perturbative computations are well known. 
LPT is hard and as a result computations of renormalization constants are usually at one loop, with 
second order computations just now making their entrance 1 . Due to the bad convergence properties 
of LPT, large use is made of so-called Boosted Perturbation Theory (BPT), often in the Tadpole- 
Improved (TI) variant. By making use of the Numerical Stochastic Perturbation Theory (NSPT) 
method 2 , we are in a position to perform quite high orders computations: this means three (and 
even four) loops for quark bilinears. By making use of BPT, we are in a position to assess con- 
vergence properties and truncation errors of the series and to make contact with non perturbative 
results, where available 3 . This is relevant, since both perturbative and non perturbative compu- 
tations rely on the same assumptions, which strictly speaking are only proved in the perturbative 
framework. 



2. Computational setup 

We performed our computations for Wilson gauge action and (unimproved, at the moment) 
Wilson fermion action (r = 1). The quenched computations have been performed on lattice sizes 
of both 32 4 and 16 4 . Unquenched computations are at the moment 4 for nj = 2,3 on 32 4 . The 
bigger lattice size fits well on our APEmille crate, while standard PC's are enough for the smaller 
size. In the following we will focus on the tif = 2 results. Configurations have been stored in order 
to perform many other computations. 

We work in the RI'-MOM scheme [||]. Quark bilinears operators are computed between (off 
shell) quark states of momentum p and then amputated to get T-functions 

( P \WV\P) = G(p) -> Tip) =S-\p)G(p)S-\p). (2.1) 

By making use of a convenient projector Po one then projects on the tree level structure 

0(p)=Tr(P T(p)) (2.2) 

to obtain the operators in terms of which the renormalization conditions are given 



ZoZ^O(p) 

Z q is the quark field renormalization constants defined via 

. 1 ^pS- 1 ) 



1. (2.3) 



Z, 



12 p 2 



(2.4) 



1 See for example [Q] at this conference. 

2 The method is reviewed in |Q| and new technical details were also presented at this conference 
3 This is the case for the quantities at hand: see |^] at this conference. 
4 ny = 4 computations are also on their way. 



237/2 



Wilson fermions quark bilinears to three loops 



Francesco Di Renzo 



Renormalization conditions are given in the zero quark mass limit in order to get a mass indepen- 
dent scheme. We note with this respect that we are in a position to stay at zero mass by plugging 
in the Dirac operator the convenient (critical mass) counterterms. These are analytically known at 
one and two loops [Q], while third loop has been computed in Going to four loop, we obtained 
as a byproduct the corresponding new order for the critical mass 5 . 

The gauge is fixed to Landau, which determines the one loop quark field anomalous dimension 
to vanish. 

One should keep in mind that the RI'-MOM scheme is defined in infinite volume, while we are 
of course forced to finite volume approximations: as we will see, care is to be taken of this point. 

In a traditional perturbative computation fixing the divergences (i.e. the log's) is the easy 
part of the job, while obtaining finite parts is usually much harder. In NSPT it is just the other 
way around. We actually take the anomalous dimensions y's for granted. In RI'-MOM they are 
known to three loop order [§] : notice that is the reason why we stop at third order. Only for RG- 
invariant quantities (like the ratio Z p /Z s on which we will focus in the following) we have almost 
no conceptually severe limitations. 

In the end, at a given loop L this is what one expects for the coefficient of a renormalization 
constant: 

L 

ZL = c L + Y d d i {y)\og{p) i + F{p) (p = pa). (2.5) 

The first (finite) term is the one we are interested in (in the continuum limit), log's coefficients are 
known function of anomalous dimensions while the remainder is a scalar function of hypercubic 
invariants. The way to extract the continuum limit from hypercubic-symmetric Taylor expansions 
is explained in [|5|], which also contains a practical example, namely for the quark field renormal- 
ization constant. 



It is useful to see how our master formula Eq. ( |2.3[ ) reads in a specific case, i.e. for the scalar 
current ( first order 6 ) 

which determines 

z^=z?-(o^-^hog(p 2 )) (2.7) 



Notice that is what is actually coming from measurements, while Zq comes from an indepen- 
dent computation (there is no log because we are in Landau gauge). Higher orders simply requires 
some more algebra in order to solve for the unknown quantities (modulo some remarks we will 
make later). 

3. Z p /Z s to four loops 

The perfect quantity to compute is the ratio Z p /Z s . Quark field renormalization constants drops 

5 A quite surprising result is that resumming the critical mass at four loop in a convenient BPT scheme one obtains 
a result which is not so far from non perturbative determinations. 



6 We slightly change notation with respect to Eq. (2.3): now the order is understood as the superscript, while sub- 
script denotes the observable. 
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Figure 1: The resummation of Z p /Z s and of Z s /Z p (tif — 2) at j3 = 5.8. On the x-axis the values of the 
different couplings xq,x\ ,X2,x?, (see text). Color code for y-axis as follows: blue is one loop result, green 
two loops, red three loops and pale blue four loops. 



out in the ratio. There is no anomalous dimension around and so there is no conceptual limitation 
to our computations (apart from practical ones, at moment mainly dictated by computer memory 
limitations). As an extra bonus, in our perturbative approach Z p /Z s and Z s /Z p come from different 
signals, so that one has the handle of verifying that the series are actually inverse of each other. 
For rif = 2 we obtain 

Z p /Z s = 1 - 0.487(2) p- 1 -1.46(2) /T 2 - 5.36(7) /3~ 3 + ... 

Z s /Z p = l+0.487(2)/r 1 + 1.70(2)/T 2 + 6.9(l)/T 3 + ... (3.1) 

One can verify that series are actually inverse of each other to a very good precision and we 
point out that finite size effects have been proved to be well under control by comparing results on 
different lattice sizes. These results come from the (quite huge) collection of three loops configu- 
rations that we stored. We actually pushed the computation even to four loops, even though with 
less statistics. At four loop we only quote the resummed result. For this ratio the computation is 
anyway safe. For quantities involving anomalous dimensions things are more cumbersome: one 
has also to fit an unknown y, which comes from log's. Next step is now to resum the series. We 
now make use of BPT and resum (for ap = 1) using different expansion parameters according to 
the definitions 

xo^/T 1 *i = ^p x 2 = ~\og{P) *3 = ^ (3.2) 

where P is the plaquette 7 . We resum at j3 = 5.8 in order to make contact with the results of [^]. 
Resummation is plotted in Figure 1 : different colors are for different orders of resummation. We 

7 While X2 and X3 are widely used, x\ was defined just in order to better point out the general picture that emerges 
from the choice of different couplings: see later. 
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Figure 2: Computation of one loop renormalization constant for the scalar current. Upper points are the 
unsubtracted o s , while lower (ticker symbols) stand for the subtracted op — log (/3 2 ). Analytic result 
is marked with a darker symbol. On the left: no correction for finite volume. On the right: finite volume 
tamed-log taken into account. 



regret that there was an error in the four loop resummation shown at the conference. We here 
display the correct four loop result, which does not change the overall picture. 

From the figure one can better understand the effect of one loop BPT: one is actually sitting on 
a straight line, whose slope is dictated by the first loop coefficient. Only if higher loops are taken 
into account a reliable understanding of what is going on can be achieved. In particular, one can 
easily force the series to oscillate widely by changing the coupling. We stress that convergence 
properties should be assessed on a case by case basis, i.e. depending on the observable. The same 
holds for TI-BPT and in general if also the scale is changed while switching to a new coupling. 

As expected, at a fixed coupling we get milder and milder variations by switching on higher 
and higher orders, while at a fixed order we get milder and milder variations by changing the 
coupling. Convergence properties of the series seem reasonable for %2 and X3. For these couplings 
it is interesting to point out that if one adds to the result at a given order the deviation from the 
immediately lower order, one always ends up at the same result, which for example in the case of 
Z p /Z s is 0.77(1) 8 . Besides that, resummed series are almost inverse of each other. The comparison 
with the non perturbative result in [^] yields encouraging result: things are working pretty well in 
the RI'-MOM scheme. 



4. Treating anomalous dimensions 



Let us go back to Eq. (2.7) (the scalar current at one loop). We make use of this one loop ex- 



ample since the analytic result is known [|3p. A naive application of Eq. ( J2.7[ ) yields what is plotted 
in Figure 2 on the left: the result is missed, as if one were subtracting too much by subtracting the 
log. The effect is systematic: a similar picture is obtained for example for Z p . By computing on 
different lattice sizes, one realizes that is a finite size effect, getting worse and worse as the volume 



This is a popular way to pin down a truncation error. 
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decreases. We stress that this is a pL == ^Na = pN effect, L being expressed in terms of the 
lattice spacing as L = Na. One could say that the log is tamed by finite volume. Since we want to 
regard our finite volume computations as inifite volume approximations, we need to correct for this 
effect. A solution comes from computing the expected log in the continuum at the same values of 
pL = pN one is interested in. The result is shown in Figure 2 on the right. Notice anyway that if 
one stays away from the lower momenta results are anyway quite safe. 

5. Conclusions and prospects 

NSPT is a valuable tool to compute logarithmic-divergent renormalization constants for lattice 
QCD. Now that configurations are stored, a lot of computations are possible and, in particular, a 
full account of quark bilinears to three loops will be ready soon. The main caveat when making 
use of our method is that care is needed when dealing with log's, because of finite volume effects. 

We are in a position to assess convergence properties of the series and to shed some light 
on BPT, a blind application of which is never advisable. Even when a non perturbative result is 
available, our method can provide an important check to better assess systematic errors. 

References 

[1] G. Martinelli, C. Pittori, C.T. Sachrajda, M. Testa and A. Vladikas, A general method for non 

perturbative renormalization of lattice operators, Nucl. Phys. B 445 (1995) 81 [hep-lat/941 1010]. 

[2] M. Luescher, R. Narayanan, P. Weisz and U. Wolff, The Schroedinger functional: a renormalizable 
probe for non abelian gauge theories, Nucl. Phys. B 384 (1992) 168 [hep-lat/9207009]. 

[3] Q.Mason, High precision fundamental constants using lattice Perturbation Theory, these proceedings, 
PoS(Lat2005). 

[4] F. Di Renzo and L. Scorzato, Numerical StochasticPerturbation Theory for full QCD, JHEP 04 (2004) 
073 [hep-lat/0410010]. 

[5] V. Miccio, F. Di Renzo, A, Mantovi, C. Torrero and L. Scorzato, Fermionic observables in Numerical 
Stochastic Perturbation Theory, these proceedings, PoS(Lat2005)108. 

[6] D. Becirevic, B. Blossier, Ph. Boucaud, V. Gimenez, V. Lubicz, F. Mescia, S. Simula and 

C. Tarantino, Non perturbatively renormalized light quark masses with two dynamical fermions, these 
proceedings, PoS(Lat2005)079 [hep-lat/0509091]. 

[7] E. Follana and H. Panagopoulos, The critical mass of Wilson fermions: a comparison of perturbative 
and Monte Carlo results, Phys. Rev. D 63 (2001) 017501 [hep-lat/0006001]; 

S. Caracciolo, A. Pelissetto and A. Rago, Two-loop critical mass for Wilson fermions, Phys. Rev. D 64 
(2001) 094506 [hep-lat/0106013]. 

[8] J. A. Gracey, Three loop anomalous dimension of non- singlet quark currents in the RF scheme, Nucl. 
Phys. B 662 (2003) 247 [hep-ph/03041 13]. 

[9] G. Martinelli and Yi-Cheng Zhang, The connection between local operators on the lattice and in the 
continuum and its relation to meson decay constants, Phys. Lett. B 123 (1983) 433. 



237/6 



